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Abstract 

We present a dispersive method which allows to investigate the low-energy couplings 
of chiral perturbation theory at the next-to-leading order (NLO) in the 1/Nc expansion, 
keeping full control of their renormalization scale dependence. Using the resonance chiral 
theory Lagrangian, we perform a NLO calculation of the scalar and pseudoscalar two- 
point functions, within the single-resonance approximation. Imposing the correct QCD 
short-distance constraints, one determines their difference Il{t) = ns(i) — Ilp{t) in terms 
of the pion decay constant and resonance masses. Its low momentum expansion fixes then 
the low-energy chiral couplings Lg and C38. At ^0 = 0.77 GeV, we obtain Lg(/Uo)'^^^^'' = 
(0.6 ± 0.4) • 10-3 and C^^sit^of^^^^ = (2 ± 6) • 10"^ 



1 The large— Nc limit 



In recent years we have witnessed a spectacular progress in our understanding of low-energy 
effective field theories [H [21 [3] . In particular, chiral perturbation theory (xPT) has been estab- 
lished as a very powerful tool to incorporate the chiral symmetry constraints when analysing the 
strong interactions in the non-perturbative regime [H [5[ E] . The precision required in present 
phenomenological applications makes necessary to include corrections of 0{p^). While many 
two-loop xPT calculations have been already performed [3 El |9l llUl 111], |12l |13l [Til [15] , the large 
number of unknown low-energy couplings (LECs) appearing at this order puts a clear limit to 
the achievable accuracy |T6] . 

The limit of an infinite number of quark colours has proved very useful to bridge the gap 
between xPT and the underlying QCD dynamics [T71 |TH]. Assuming confinement, the strong 
dynamics at Nq — > cxd is given by tree diagrams with infinite sums of hadron exchanges, which 
correspond to the tree approximation of some local effective Lagrangian fi9[ [20] . Resonance 
chiral theory (RxT) provides a correct framework to incorporate the massive mesonic states 
within an effective Lagrangian formalism [21]. Integrating out the heavy fields one recovers 
the xPT Lagrangian with explicit values of the chiral LECs in terms of resonance parameters. 
Moreover, the short-distance properties of QCD impose stringent constraints on the low-energy 
parameters [22] . 

Truncating the infinite tower of meson resonances to the lowest states with ^, 1 
and l"*"^ quantum numbers (single-resonance approximation, SRA), one gets a very successful 
prediction of the 0{p^Nc) xPT couplings in terms of only three parameters: My, Ms and 
the pion decay constant F [T7]. Some O(p^) LECs have been already predicted in this way, 
by studying an appropriate set of three-point functions [23]. More recently, the program to 
determine all 0{p^) LECs at leading order in 1/Nc has been put on very solid grounds, with a 
complete classification of the needed terms in the R^T Lagrangian |24j . 

Since chiral loop corrections are of next-to- leading order (NLO) in the 1/Nc expansion, the 
la.Tge-Nc determination of the LECs is unable to control their renormalization-scale dependence. 
For couplings related with the scalar sector this introduces large uncertainties, because their fi 
dependence is very sizable. A first analysis of resonance loop contributions to the running of 
L^Q(/i) was attempted in Ref. [25]. More recently, a NLO determination of the xPT coupling 
Lg{n) has been achieved, through a one- loop calculation of the vector form factor in R^T |26j . 
In spite of all the complexity associated with the still not so well understood renormalization 
of RxT [2B| [27] . this calculation has shown the potential predictability at the NLO in 1/Nc- 

In this article we present a NLO determination of the couplings Ll{^) and Cl^{^). Using 
analyticity and unitarity we can avoid all technicalities associated with the renormalization 
procedure, reducing the calculation to tree-level diagrams plus dispersion relations. This allows 
to understand the underlying physics in a much more transparent way. In particular, the subtle 
cancellations among many unknown renormalized couplings found in Ref. [2B] and the relative 
simplicity of the final result can be better understood in terms of the imposed short-distance 
constraints. 

Let us consider the two-point correlation functions of two scalar or pseudoscalar currents, 
in the chiral limit. Of particular interest is their difference n(t) = n5'(t) — np(t), which is 
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identically zero in QCD perturbation theory. When t oo, this correlator vanishes as 
with a coefficient proportional to as{qTqqTq) [25112^ . The low-momentum expansion of n(t) is 
determined by xPT to have the form [U [7] 

f r 

m = 5oM ^ + 32L^(/.) + 



t 

+ ^ 



32C3^3(/.)-I^h-ln4)+O(iV0)J +0(t^) [, (1) 




with Tg = 5/48 [3/16] and vf^ = -5Lg/6 [-3ivi:/2] in the SU{3)l ® 5^/(3)r [U{3)l ® ?7(3)r] 
effective theory [HI?]. The correlator is proportional to Bq = {qq)^/F'^, which guarantees the 
right dependence with the QCD renormalization scale. The couplings F^, Lg and C^s/F'^ are 
0{Nc), while Tg and T^^^/F^ are 0(1) and represent a NLO effect^ 
In the large- A^^c limit, n(t) has the general form 



^2 8rf2. 



which involves an infinite number of scalar and pseudoscalar meson exchanges. This expression 
can be easily obtained within R^T, with Crm and drm being the relevant meson couplings. For a 
finite number of resonancesjl one finds that imposing the right high-energy behaviour (~ 1/t^) 
constrains the resonance parameters to obey the relations [33] : 

E (cm, - dl^) = ^ , E {^Ih Ml - Ml) =5, (3) 

i i 

where 6 = SttOsF'^/A ^ 0.08 OgF^ x (1 GeV)^. Truncating the infinite sums to their ffist con- 
tributing states and neglecting 6, these relations fix the corresponding scalar and pseudoscalar 
couplings in terms of the resonance masses: 

. = Mp = Ms 

8 M2-M|' 8 M2-M|' ^ ^ 

Note that Eqs. (j4]) imposes Mp > Ms- On the other hand, the low-energy expansion of ([2]) 
determines [22| [24] 

+ 
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Using the approximate constraint Mp/y/2 ^ Ms ~ 1 GeV [34j, this gives Lg ^ 3F^/(32M|) ^ 
0.7 ■ 10^^ and Csg ~ 7F^/(64M|) ^ 7 ■ 10^^. However, one does not known at which scale /x 
these predictions apply. 



^Note that we have defined a dimensionless C^s and C^g^; in the notation of Ref. [7] that corresponds to F^G^s. 
and respectively. 

^ Some issues related to the truncation of the spectrum to a finite number of resonances are discussed in 
Refs. |30l[3Tl[32]. 
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2 NLO corrections 



At the NLO in l/Nc, n(t) has a contribution from one-particle exchanges, with the structure in 
Eq. ([2]), plus one-loop corrections An(t) generating absorptive contributions from two-particle 
exchanges. The corresponding spectral functions of the scalar and pseudoscalar correlators take 
the form: 



with 



-imn5(t) 

TT 



TT 



Imnp(t) 



8c^5(t-M|) 



167r2 J ' 



F^5(t)+8ci^5(t-M|) + ^| 



(6) 



Ps{t) 



ppit) 




\Fi\t)\' 



4M^ 

eit - Ml) (i - 



1^1" wr +■■ 



(7) 



(8) 



We have adopted the single-resonance approximation and, moreover, we have only taken explic- 
itly into account the lowest-mass two-particle exchanges: two Goldstone bosons or one Goldstone 
and one heavy resonance. In the energy region we are interested in, exchanges of two heavy res- 
onances or higher multiplicity states are kinematically suppressed. Our normalization takes into 
account the different flavour-exchange possibilities. The relevant two-particle cuts are governed 
by the following scalar. 
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and pseudoscalar. 
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form factors |35j . 

The RxT couplings [211 [23] generating these form factors have been determined imposing a 
good high-energy behaviour of the corresponding spectral functions, i.e. that the individual form 
factor contributions to ps{t) and pp{t) should vanish at infinite momentum transfer. Moreover, 
we have used the constraints (jl]) and the analogous relations (Weinberg sum rules and good 
high-energy behaviour of the vector form factor) emerging in the vector /axial- vector sector. It 
is quite remarkable that these short- distance constraints completely determine the form factors 
in terms of the resonance masses [35j. The form factor F^'^{t) turns out to be identically zero, 
within the SRAH 

Using its known analyticity properties, AIl(t) can be obtained from the spectral functions 
through a dispersion relation, up to a subtraction term which has the same structure as the 
tree-level scalar and pseudoscalar resonance exchanges. Therefore, the unknown subtraction 
constants can be absorbed by a redefinition of Cm, dm, Ms and Mp at NLO in l/Nc'- 

. + — + An(t)L \ . (11) 

The explicit expression of An(t)|p is relegated to appendix A. At large values of t, the one-loop 
contribution has the behaviour 

F2 F^Mi f^,,, , -t\ „/l 

Since the logarithmic term ln(— t)/t^ should vanish, one obtains the constraint 



e^:^l«fl-^l^ + 3fl-2t)No, ,13) 



leading to 



327r2F2 \ \^ Ml) Ml \ Af|^ 



M2\ M| / M| 



MlJ Ml\ 2Ml^ 



(14) 



which requires Ma < \/2My. Imposing the right short-distance behaviour (~ 1/t^) in n(t), one 
gets 

F"Ml5^^l - Scl^'M^s' + ^'CM^p' = -86, (15) 

where the corrections 

Vo - 32^2^2 Mj,)^'^ ^^[mI V Ml y- Ml) / ^^^^ 



^We are using the SRA and only operators constructed with up to 0{p^) chiral tensors are aUowed in the 
R^T Lagrangian. While the inclusion of extra resonances would certainly generate a non-zero value for F^'^{t), 
its numerical impact in our final resuls would be small, because the short-distance contraints would force a 
compensating effect with the other contributions. The same comment applies to possible higher-derivative chiral 
structures, which generate a bad high-energy behaviour and, therefore, would need to be severely tuned. 
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are known functions of the resonance masses 




Note that Eqs. (ITS!) determine the NLO couphngs c 
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3 L^(/i) at NLO 



The low-momentum expansion of the R^T correlator in Eq. (ITTi) reproduces the U{3)l (S> U{3)r 
xPT result ([1]), with a definite prediction for the LEG LKfi): 

,2 



with 
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(19) 



(20) 



We have used the relations in Eqs. (fTSl) to ehminate the exphcit dependence on the effective 
couphngs and d^^. 

Eq. (1191) modifies the large- iVj^ result in with NLO corrections S^^^, S^^^ and A, which are 
fully known in terms of resonance masses. We have also taken into account the tiny correction 
S. Moreover, our calculation has generated the right renormalization-scale dependence, giving 
rise to an absolute prediction for the scale-independent parameter Lq^^\ Since we are working 
within the large- A^^c framework, the Goldstone-nonet loops reproduce the non- analytic In (— t) 
structure that arises in U{3)l Cg) U{3)ji xPT. To make contact with the usual SU (3) l ® SU {3)r 
theory, we still need to integrate out the singlet rji field. Computing the massive one-loop rji 
contribution to Il{t), one easily gets the known relation [36] between the corresponding Lg 
couplings in the two chiral effective theories. At O(p^) in the U{3)l f/(3)ij theory, one finds: 

The different input parameters are defined in the chiral limit. We take the ranges [H [T71 
EZl ESI ES] My = (770 ± 5) MeV, = (1.14 ± 0.16) GeV, = (1.3 ± 0.1) GeV, M^^ = 
(0.85 ± 0.05) GeV and F = (89 ± 2) MeV, and use Eq. to fix Ma, keeping the constraint 
Mp > Ms from Eqs. (jl]) and imposing Ma > 1 GeV. The correction 6 turns out to be negligible. 
One obtains the numerical prediction 

Lg^^^^ = (0.4 ± 0.4) ■ 10"^ . (22) 

The largest uncertainties originate in the badly known values of Mg and Mp, which already 
appear in the leading order prediction ([5]). To account for the higher-mass intermediate states 
which have been neglected in ([2]), we have added an additional truncation error equal to 0.12 ■ 
10~^, the contribution of the heaviest included channel (-Pvr). All errors have been added in 
quadrature. At the usual xPT renormalization scale fiQ = 0.77 GeV, Eq. (!22l) implies 

Llil^of^^^^ = (0.6 ± 0.4) ■ 10-3 , (23) 

to be compared with the value Ll(fio)^^^^^ = (0.9 ± 0.3) • 10"^, usually adopted in 0{p^) 
phenomenological analyses, or LKhq)^^^^^ = (0.62 ± 0.20) • 10"'^, obtained from the 0{p^) fit of 
Ref. dl]. 

The sizable numerical difference between L^'^^^^ and LKhq)^^^^^ shows the large sensitivity 
of this coupling to the xPT renormalization scale. This is a general trend for those LECs which 
are dominated by scalar or pseudoscalar resonance exchanges. Therefore, to perform accurate 
phenomenological applications one needs to control the renormalization scale dependence, which 
requires a determination of the xPT couplings at NLO in 1/Nc, like the one presented here for 
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4 C^8(/i) at NLO 

Following the same procedure used in Section 3, one is able to find a NLO prediction of the 
LEG CM-- 
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(25) 



We have used again Eqs. (fT8|) to fix cj!„ and (i^. As in the case of Ll{^), one has to make contact 
with the usual SU{3)l ® SU{3)r theory by computing the massive one-loop rji contribution to 
n(t). It is straightforward to get the expression that relates C^s in both theories: 
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(26) 



where we have used the LO prediction of the C(p^) chiral coupling L5 = CdCm/Ml 

F'^/{ml) [21] • 

Taking the same input parameters than in the previous section, one gets the numerical 
prediction 



C- 



5(7(3) 
38 



-1 ± 6) ■ 10" 



(27) 



At the usual xPT renormalization scale /iq = 0.77 GeV, Eq. (p7|) gives 

Q8(/^o)^^^'^ = (2±6)■10-^ 
showing again the large numerical sensitivity to the choice of scale. 
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A sizable difference between the phenomeno logical value of some 0{p^) chiral couplings and 
their large- iVc estimates was pointed out in Ref. [40|. Our result shows how the NLO corrections 
in I/Nq become relevant in some cases. These NLO contributions must be considered for a 
proper determination of the LECs. Our calculation reproduces the correct scale-dependence of 
Cl^[li) at the NLO in I/Nq. However, the 0{p^) LECs contain additional dependences on /i, 
which are suppressed by two powers of 1/ Nc [Z] . In order to remove the corresponding numerical 
uncertainty, it would be necessary to perform the calculation at the next-to-next-to-leading order 
in the 1/Nc expansion. 

5 Conclusions 

The large-iVc limit provides a solid theoretical framework to understand the success of resonance 
saturation in low-energy phenomenology [I7j. However, this limit is unable to pin down the 
scale dependence of the xPT couplings. Although this is a NLO effect in the I/Nq expansion, 
its numerical impact is very sizable in couplings which are dominated by scalar or pseudoscalar 
exchanges. 

In this paper we have presented a NLO prediction of the 0{p^) coupling L^fi), which exactly 
reproduces its right renormalization-scale dependence. Moreover, we have also determined the 
0{p^) coupling C3g(/i) at the NLO, controlling its /i dependence up to small NNLO effects. 

The determination of this two LECs has been made possible through a NLO calculation 
in 1/Nc of the H(t) = H5(t) — Hp(t) correlator, in the chiral limit. We have used the R^T 
Lagrangian, within the SRA, to compute the one and two-particle exchange contributions to 
the absorptive part of the correlator. It is remarkable that, imposing a good short-distance 
behaviour for the corresponding scalar and pseudoscalar spectral functions, one fully determines 
the relevant contributing form factors. Using a dispersion relation, we have reconstructed the 
correlator, up to a subtraction term which has the same structure as the tree-level one-particle 
contributions. 

The stringent short-distance QCD constraints on H(t) have allowed us to fix the subtraction 
constants in terms of resonance masses, with the results shown in Eqs. f lTSj) . Therefore, we have 
obtained a complete NLO determination of H(t), which only depends on the pion decay constant 
F and the meson masses My, Ma, Ms and Mp. Its low momentum expansion reproduces the 
right xPT expression, with explicit values for the LECs L^fi) and C^glfi). 

Integrating out the rji field, one can further connect the U{3)L^U{3)ji and SU (3) SU (3) ji 
effective theories. This introduces an additional dependence on M^^. At the usual scale /iq = 
0.77 GeV, we finally obtain the numerical predictions: 



The ideas discussed in this article can be applied to generic Green functions, which opens 
a way to investigate other chiral LECs at NLO in the laige-Nc expansion. Further work along 
these lines is in progress [35j . 




Ll{f^of''^''> = (0.6±0.4)■10-^ 



C^si^^of""^'^ = (2±6).10-^ 



(29) 
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A The one- loop correction AIl{t)\p 

In this appendix we show the explicit expression of the one-loop correction AIl(t)\p, generated 
by the considered two-particles exchanges, which has been calculated by using the dispersive 
method discussed in Section 2. 



An,_,(t)|v. 



An,_,(t)|p^ 
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